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O ■ We study radion stabilization in the compact Randall- Sundrum model by intro- 
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—^ I taking into account the backreactions from the scalar field on the metric. Our gen- 

QQ ■ eralization reconciles the radion potential found by Goldberger and Wise with the 
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C^ ' in the regime where the 1/A^ expansion is justified, the gravitational wave signal is 
detectable by LISA. 
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I. INTRODUCTION 

Even though the Standard Model (SM) of particle physics has achieved many impressive 
experimental successes it fails to provide an explanation to some experimental and theoretical 
issues. On the experimental side, neither the observed dark matter density can be explained 
(the SM lacking an appropriate natural candidate for it) nor the baryon asymmetry of the 
Universe can be accounted of, mainly because the SM does not provide sizeable CP violation 
sources and a strong enough first-order electroweak phase transition. On the theoretical 
side there is no plausible explanation for the huge hierarchy between the electroweak scale, 
responsible for the mass of the weak gauge bosons, and the Planck scale, apparent in the 
weakness of gravitational interactions. Thus an extension of the SM seems necessary. 

Guided by the naturalness criterion to solve the Higgs hierarchy problem one very at- 
tractive possibility is Randall-Sundrum (RS) models, which are based on the framework of 
a compact warped extra dimension [l] with two branes localized on it: an ultraviolet (UV) 
brane, which provides the UV cutoff of the higher dimensional theory, and an infrared (IR) 
one, which spontaneously breaks the conformal symmetry of the theory. In this class of mod- 
els the four-dimensional part of the metric has a strong dependence on the fifth coordinate, 
which is unobservable macroscopically, and the hierarchy between the two energy scales is 
generated by localizing the relevant physical degrees of freedom responsible for electroweak 
breaking at (or near) the IR brane. In this way the huge hierarchy can be explained by a 
natural distance between the two branes of (9(10) times the fundamental five-dimensional 
Planck length. 

To complete the picture the distance between the branes should not be considered as a 
fundamental input but instead it should arise from a stabilization mechanism. In fact this 
is even essential to avoid a massless radion generating an unobserved fifth force, and to 
obtain the observed Friedman- Robert son- Waker cosmology at late times 2|, |3| . One elegant 
possibility to stabilize the brane distance is to assume the existence of a bulk scalar with a 



five-dimensional mass ^ that is slig 
pioneered by Goldberger and Wise 



itly smaller than the fundamental Planck mass scale, as 
6|. In this setup one usually assumes the weak field limit 



in which case the dynamics of the bulk scalar field and the metric decouple. This does not 



^ Alternatively, Casimir effect can also lead to a stabilizing potential [J, l5|. 



only have the advantage of siniphfying the analysis but also facilitates the interpretation of 
the scalar action as a radion potential. 

In the modern context of the AdS/CFT correspondence [TJ] this kind of model can be 



interpreted as dua^ 
a UV completion 



to a strongly coupled gaioge theory that for instance might serve as 
8| to little Higgs models [9|. Usually the five-dimensional system is 
considered in the limit of a large ratio between the Planck mass and the bulk cosmological 
constant, which implies a large number of degrees of freedom in the dual field theory and 
allows to neglect stringy effects. However there is a certain tension between this assumption 
and a viable cosmology: even though the model has the above discussed brane setup at 
low temperature (called RS-GW in the following), at high temperature it is represented 
by an AdS-Schwarzschild (AdS-S) bulk metric according to the AdS/CFT correspondence. 
Even though a first-order phase transition could allow the Universe to escape from the 
(conformal) AdS-S to today's (quasi-conformal) RS-GW phase [lO|, its completion leads to 
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Il4l | which jeopardizes the consistency of the 



a stringent constraint on the above ratio 
original assumptions. 

On the other hand several appealing features characterize this necessary phase transition. 
One of them is its extreme supercooling, which may couple the conformal phase transition 
to the electroweak one making the latter strong enough to plausibly explain the baryon 
asymmetry of the Universe 13|. Besides a further feature is its potential testability: the 
phase transition could mark the gravitational wave spectrum which could allow to probe 
most of the interesting parameter space of the model |ll|. 

In conclusion given the success of RS models as extensions of the SM we believe it is 
of importance to clarify the impact of the phase transition on the validity of the models. 
This paper is dedicated to this issue in the case of RS-GW models. Our aim is to alleviate 
the tension between the ratio of the Planck mass to the bulk cosmological constant and the 
phase transition completion and to show how this scenario can be tested in the gravitational 
wave spectrum. 

In this paper we will proceed as follows. After introducing some notations and conventions 



in sec. |TT] we will review in sec. H 



T] the usual procedure to determine the effect of the bulk 



scalar on the brane distance [6|, |15| . We will remark that the method of ref . [6| determines the 
effective potential of the radion although without taking into account the backreactions of 
the scalar field on the gravitational metric. On the contrary the superpotential method 



solves the problem exactly but it does not provide the radion effective potential which is 
needed to study the phase transition. For this reason in sec. II VI we will present an alternative 
approach to determine the radion potential in the regime of detuned brane tensions and 
sizeable backreactions. It is based on fitting the radion potential information that we can 
determine in this regime: the position of the radion potential extrema and the radion mass 
and cosmological constant in these extrema. The obtained potential will be used in sec. |V] 
to study the phase transition and we will find relaxed bounds that may alleviate the above 
problematic parameter tension. The corresponding gravitational wave spectrum will be 
determined in sec. I VI I and we will conclude by commenting on the prospects of detection 
with forthcoming experiments as LISA. Finally we will devote sec. IVIII to summarize the 
main results of the paper and we leave in appendix |A] some technical details that we need 
in sec. llVl to determine the radion mass in the presence of a cosmological constant. 

II. NOTATION AND CONVENTIONS 

A very interesting feature of RS models with a scalar in the bulk is that the scalar field 
can stabilize the brane distance p, |l5|. The corresponding five-dimensional (5D) action is 
given by 

S= fd'x^\detgMN\ -M'R+^{d<pf-V{<f)) - $^ / ci^x^l det^^,|A„(0) , (1) 

where M is the 5D Planck scale, Aq, and V are the brane and bulk potentials of the scalar 
field 0, and the metric qmn is defined in proper coordinates by 

ds^ = e^^^'^gf^^dx^dx" - dr^ , 

Mi : g^^ = ri^^ , (2) 

dSi : g^udx^dx"" = dt^ - e^'^\dx\ + dxl + dxl) , (3) 

AdSi ■ g^Jx^'dx" = -dxl - e'^'^^''' {dxj + dxl - dt^) ■ (4) 

Accordingly the induced four-dimensional (4D) metric cjf^,^ is Minkowski, de Sitter or anti-de 
Sitter and the corresponding 4D cosmological constant A has mass dimension equal to 2. In 



all cases the Ricci-tensor and scalar turn out to be 



R,, = e^\AA' - A' - 3Ae'^)-g,, , (5) 

i?55 = -4A'2 - AA" , (6) 

R = 20A'^ + 8A" - 12Ae-2^ , (7) 

where ' = d/dr. The equation of motion for the scalar field reads 

a 

and the Einstein equations have the form 

A" + Ae-'^ = -^r - y E ^"^(^ " ^-) ' (9) 

a 
^2 _ ^-2A = _'^V+ — 0'2 , (10) 

with K^ = 1/(2M^). Hereby the localized terms impose the following constraints (assuming 
a Z2 symmetry across the branes) 

A l^_^ = -yA„(0(r„)) , |,^_^ = — . (11) 

Using these equations in the action one obtains 



S = QM-' / rf'xVl det c/Afjv|e"'^A . (12) 

We would like to emphasize that this relation does not rely on any approximation so far. 
The key observation is that from an effective four- dimensional point of view the expansion 
parameter A has to be related to the value of the radion potential. If the system fl8|)- flT0|) 
allows for several solutions, the corresponding values of A allow to determine the difference 
in potential energy of the radion between different configurations. Here we will focus on 
positive cosmological constants and our aim will be to determine A with an accuracy that 
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goes beyond the usual weak field assumption cf? <^ M^ |6| 



III. STABILIZATION AND GW MECHANISM 

Solving the system (!8|)-( IT0|) is for generic scalar potentials a hard task. A possibility to 
overcome this difficulty is to consider the quadratic Goldberger-Wise (GW) scalar poten- 



tial [61 

^(0) = -^ + ^-^^ , (13) 

Aa(0) = A°+7„(0-t;„)2 with 7„ ^ oo , (14) 

where / is the AdS length such that 1// = A; is of the order the Planck scale, and v^ have 
mass dimension 3/2 and 7q, has mass dimension 1. We assume that these potentials are 
chosen such that the warping in the metric is close to AdS, meaning that A{r) in flTOj) is 
dominated by the bulk vacuum energy 

/20'2 < 24M^ fm^cj)'^ < 24M^ A < T^ . (15) 

At leading order in the scalar sector then decouples completely from gravity. The equations 
of motion for the metric are then given by A = and A = —1// and the gauge choice 
A{0) = leads to the solution 

A{r) = -r/l , (16) 

in the bulk. We assume the UV and IR branes to be located at ri = and r2 = Tq 
respectively. 

Finally under the assumptions flTSj) the solution of the scalar field turns out to be 

0(r) = v,e'- + (v, - v,e'-^) fZ~S-r. ' ^^^^ 



lk± = 2± V4 + m2/2 , (18) 

where Va denote the values of at the UV and IR branes ^. We focus on negative values ^ 
of /c_ and a solution of the hierarchy problem (ro ~ 37/) requires 

t;i ~ f 2 ~ M3/2 , 1 ~ e'=-''« < 10^^ ~ e''"/' , (19) 

such that one obtains 

0(r) ^ i;ie'=-" + {v2 - y^e'=-"«)e'=+("-"o) . (20) 

Using this in the action ([T]) yields 

S = r^ I d^x hk.vl - e"^"°/' [(4 - lk.){v2 - vie''-''f + lk.vi\ } . (21) 



^ In the brane potentials (ITil) the hnear terms satisfying eqs. PT|) and ITTl) are omitted since the results of 

our analysis will be independent of them. 
^ Most results in our analysis are easily carried over to positive fc_ . 
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In ref. [6| this is identified (after a cliange of sign) with the potential of the radion Vqw ^^ 

S = -jd'xVGw{ro), (22) 

and has a minimum at ^ = ^_ (^.e. Tq = r_) where 

^^^gfc.ro^ ^_^!^e^-'^- = 1 + |/A:_ 1^/72- /A;_/4 + 0((/fc_)3/2) . (23) 

The value of the potential at this minimum is of (9(TeV^) 

Vcwir-) - Vgw{oo) ^ r^ \lk_f\le-^'~'' . (24) 

However this result is not conclusive. In an effective description of the radion one would 
like to split the dynamics into four-dimensional gravity and the radion degree of freedom. 
Hence the radion potential should not only depend on the scalar part of the action but 
might also receive a contribution from 5D gravity. In the derivation of eq. fl2^ we neglected 
contributions to the field A of order 0^ but these terms can potentially change the action 
and thus the potential seen by the radion. In fact from eq. ( !T2|) it is intuitive that these 
corrections should arise since, in a non-expanding background, A = 0, the action should 
vanish. In particular these additional contributions will change the boundary conditions on 
the branes and could potentially modify the difference in action between the minimum and 
the limit tq — )■ 00 which is most important for the analysis of the phase transition. 

Alternatively the system (l8l)-( TT0|) for A = can be solved by the so-called superpotential 



method 
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16| . Its large advantage is that it provides exact solutions since backreactions 
are automatically taken into account. Starting from a superpotential W^(0) the equations 
of motion for the choice 

eaA„(0) = W^(0(r„)) + ^^(^(^")) (0 _ 0(^^)) + ^^^^(0 _ 0(^^))2 ^ (25) 

[where £12 = ±1 refer to the two branes, at ri = and r2 = Tq respectively, according 
to the Z2-orbifold boundary conditions] can be recast in term of the first-order differential 
equations 



In particular the superpotential of the form 

12M3 



W = ^^ + A;_02 , (28) 



provides the following scalar potential 



1^^ 2/2 '^ 6P 



^ = J^ + ^72 ^ TH^^^ ■ (29) 



For the aim of estimating A at order (p"^, eq. f l29|) coincides with the potential f lT3|) after 
imposing eq. (|T8l) an d requiring 0^ ^ M^. The corresponding solutions to the equations of 



motion read 
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16| 



A = -J- Ue''-"- , (30) 

/ o 

(j) = Vie'"-' . (31) 

Consequently once one chooses V2 the brane distance is fixed by 

e_ = ^e'^-- = 1 , (32) 

which differs from eq. ( l23l) . Furthermore notice that this solution is based on the requirement 
A = 0, which means that the action at its extremum C, = C,~ vanishes [cf. eq. (IT2l)]. On 
the other hand in the limit tq — t- oo the (p profile f l3T]) . which matches with f lT7|) . is still a 
solution of the equations of motion with A = and thus the radion potential also approaches 
asymptotically a vanishing cosmological constant. Hence the superpotential method shows 
that backreactions can have an important impact on the radion potential. 

IV. RADION POTENTIAL INCLUDING BACKREACTIONS 

The superpotential method can also be generalized to non-vanishing cosmological con- 



stant [16| and in principle every solution to the equations (l8l)-(TT0l) can be exactly derived 
from some superpotential. Nevertheless it is not of much use in determining the radion 
potential since, for fixed scalar bulk and brane potentials, it cannot be used to find several 
solutions corresponding to different brane separations ^. Understanding the structure of the 
radion potential (partially) including backreactions is the aim of this section. 

^ Using the superpotential method a change in the brane separation would lead, according to (I25[) . to a 
change in the brane potentials. 
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For an arbitrary brane separation the system ([8])- ([9]) does not always have a solution. 
There are three integration constants, the brane separation tq and the parameter A and four 
constraints on the branes ( iTTl) . However only the combination Ae^^ enters in the equations 
and one integration constant can be eliminated. It can be used to choose e.g. A{Q) = or 
|A| = 1. In summary we have four constants to be fixed by four boundary conditions, so 
generically one expects a unique solution for given bulk and brane potentials. This is not 
too surprising since one would expect that the system does not allow for a time-independent 
solution (up to Hubble expansion) for the radion field when it is not located at an extremum 
of the potential. 

Naively one would like to determine the action for several brane separations and identify 
it with the negative potential seen by the radion. This is basically the procedure followed by 
Goldberger and Wise in a fixed gravitational background. However there are two objections 
to calculating the radion potential in this way if backreactions are included. First, solutions 
to the Einstein equations do not constitute an extremum of the Einstein-Hilbert action, 
so the action of the gravitational part should have no physical significance. This problem 
can be easily overcome by including the Gibbons-Hawking term in the action [171, [l8| as we 
will see. Second, the system of equations fl8|)- flT0|) only allows for brane separations that 
correspond to extrema in the radion potential. One way of avoiding this latter problem 
would be to solve all equations includinga time-dependence and this program is followed in 
the vicinity of the static solution in ref. [3|, though without any stabilization mechanism for 
the radion. 

In the present work we will present an alternative and simplified treatment which includes 
the bulk potential. We first determine the action and the radion mass in the extrema of 
the radion potential and, after computing the kinetic term, we use this information to get 
a reliable fit to the whole potential. We finally check that the results of our effective four- 
dimensional theory are consistent with the ones of ref. 3| . 

We will focus on the particular case of the bulk and brane potentials (fT3l) - (iT4l) constrained 
by the bounds ( IT5l) and with solution ( l20l) . However the rationale we follow can be used to 
analyze the radion potential also in other scenarios. 
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A. Cosmological constant and position of the extrema 

The values of the action at its extrema are related to the expansion parameter A at those 
points according to ( IT2l) . We will determine this parameter in this section. A basic ingredient 
for our analysis is the scalar solution ( l20l) which was obtained under the constraints (TTSll . 
Besides we make use of stiff scalar potentials on the branes. In this case the scalar field 
is fixed to the values Va on the branes and eq. fITOl) in combination with the boundary 
conditions (TTT]) reads 



12^ 






^K{va) + -V{v^) . (33) 



On the basis of this equation it is possible to accurately determine the expansion parameter A 
including ^ backreactions. Notice that due to the presence of stiff brane potentials only the 
left-hand side depends on the brane separation, once the model parameters are fixed, while 
the right-hand side can be arbitrarily chosen due to the brane potentials. In the following 
we will first discuss the case where the brane potential is tuned to reproduce the results 
obtained with the superpotential method for the choice 

,n 12M3 ,2 
SaX'a = —j- + k-vl , (34) 

and we will subsequently detune it to obtain cosmologically more realistic potentials as the 
one deduced by Goldberger and Wise. 

1. Tuned case 

Let us consider eq. (l33ll evaluated for the two branes at r G {0,ro}. For the brane 
potentials A° used in the superpotential method these two equations read 

24M=^Ai = vlk^_ - {vik. + k+{v2 - vie^-''°)e-^+'^f , 
24M=^A2e2^'« = vlkl - {vik_e''-'° + k+{v2 - Vie^-'''')f , (35) 

where we have made use of the equality 

m^ = k'i-Ak_/l . (36) 



Subject to the constraints in (IT5|) . 
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Owing to the choice of brane potentials two solutions with vanishing expansion parameter, 
A = 0, are given by tq — t- oo and tq = r_ with 



■^1 k 

V2 



(37) 



as obtained by the superpotential method. Nevertheless there is an additional solution 
with a larger brane separation and a positive cosmological constant that leads to a positive 
Einstein-Hilbert action according to eq. ( TT2l) . As an example we show in fig. [T]the functions 
Ai 2 versus vq corresponding to a given set of parameters. 

An analytic estimate for the second solution can be obtained as follows. Notice that the 
function Ai is small compared to A2 due to the factor e~^+^° and 



24M^Ai ^ -2k^k+vi{v2 - v^e 



k^ro\ -k+TQ 



(38) 



Hence the values of tq for the solutions to Ai = A2 are close to the zeros of A2 that are given 
by 



i. 



^1 fc- 



fc+ =F fc_ 



f 2 kj^ — A;_ 

The first solution is just the usual solution obtained by the superpotential method 

r_ = — log — , 



(39) 



(40) 
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FIG. 1: The functions Ai^2 leading to the two different solutions of the equations of motion. The 
parameters used are k^l = 0.5, vi = O.IM^'^, V2 = O.OSAf^'^ and/ = M~'^ which corresponds to 
m^t^ = 2.25. The values are chosen for illustrative purposes and do not lead to a realistic hierarchy. 
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while the other solution is given by 



r+ = r_ + — log- ^r_ + — ^r_ + -. (41) 



Using this value in eq. (!38l) gives for A the value 



The expansion parameter A+ is in this case additionally suppressed compared to the energy 
scale involved in the hierarchy problem and one might attempt to use this fact to solve, or 
at least to alleviate, the cosmological constant problem in a similar way to the proposals in 
ref. [l9|. However and not unexpectedly, as we will find in a later section, this solution does 
correspond to a maximum in the radion potential and is not stable. 

Let us compare this result with the findings from the superpotential method. Since back- 
reactions are taken into account to order 0^ the action and the brane separation r_ should 
agree to this order and they indeed precisely do. Notice that the expansion parameter A 
for an infinitely large brane separation also vanishes in agreement with the superpoten- 
tial method. One can thus infer that the radion potential has two degenerate minima at 
tq = {r_,oo} separated by a maximum at r+ ^. Consequently a realistic phase transition 
proceeding from oo to r_ requires some detuning away from the scenario with a superpo- 
tential, e.g. by changing the brane potentials which is the next step of our analysis. 

2. Detuned case 

Using the brane potentials of the tuned case above does not allow a realistic phase tran- 
sition. Since the radion potential is degenerate at tq = r_ and tq — )■ oo the system would 
be stuck in the high temperature phase. For a realistic model we should modify these two 
configurations by some detuning of the UV and IR brane tensions. Following the same steps 
leading to eqs. ( l35l) the new boundary conditions for A' in eq. ( ITTi) lead to 

24M3Ai = vlk^_ - {vik_ + k+{v2 - t;ie'=-"«)e-'=+"«)2 + c^vf , (43) 

24M=^A2e2'^° = vlk^_ - {vik.e^-'^ + k+{v2 - vic^-''))^ + C2vl , (44) 

where we have parameterized the detuning by Caf^ = |(A^ ~ -^a ) • 



A more detailed explanation why this guess is correct will be provided in the next subsection. 
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We first enforce tliat today's observed expansion is small. If the minimum of the potential 
has a vanishing expansion parameter A, according to eq. ( H4|) it is given by 



e_ = ^e-- = ^^|V^^V+^ , (45) 

which depends on the free parameter C2. For A = eq. ( H3|) reads 

= 24M^Ai ^ 2viV2k+k. (^ - 1) 6"'=+'^ + Cit;!^ , (46) 

which together with the position of the minimum ^_ in eq. ( H5|) fixes the value of Ci as 

cit;2 = 2v,V2K\k_\{^_ - 1) €-'=+'■- . (47) 

This choice is the fine-tuning that is needed to solve the cosmological constant problem and 
it is in general present in RS-type models. 

Once one fixes ^_ (and correspondingly Ci and C2) the system ( H3|) - (H^ has a second 
solution C,+ that fulfills 

k'. + C2- (fc-e+ + ^+(1 - e+))' ^ ci4e'''+/' 

= 2-k+\k^\{^_ - 1) e-'^+'-'+^^+Z' . (48) 

V2 

As long as C2 is small the right-hand side can be neglected and the position of ^+ is, in a 
similar way to the minimum ^_, given by 

J, = ^e-- = ^l^^^ = -i^-^. (49) 

On the other hand if C2 surpasses k\ — k'^ the left-hand side is positive for ^+ > 0. In this case, 
the right-hand side has to become comparatively large what happens at r_|_ ^ "2"'"- ~ ^^- 
which implies ^+ ~ — ^^ < ^_ for A;_ < 0. Comparison shows that sizeable deviations from 
fj49|) are hence not expected for 



2f 2 V 4f i ^^2 



?-^-^ + t/S + ^. (50) 



It turns out that the parameter space that violates this bound is less interesting in the sense 
that it leads generally to rather strong phase transitions. In this way we will mostly be 
concerned with the region in parameter space where (149!) is a good approximation. 
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Plugging eq. (H7|) in (H3|) one finds the expansion parameter at the maximum ^_|_ to be 

A+ ^ ^^+1^-1 ((e- - l)e-'^-^'- + (1 - e+)e-'^-^-) . (51) 

If ^_ is not too close to unity the first contribution dominates this expression because 
r_|_ > r_. On the other hand the cosmological constant in the limit .^ — > (ro — t- oo) is given 
by this first contribution 

as it can be easily deduced from eqs. (l43ll and (l47l) . The radion potential as a function of ^ 
is hence very flat between ^ = and ^ = C,+ and then quickly drops to zero at .^ = ^_. In 
a subsequent subsection we will provide a parametrization that interpolates between these 
extremal values. 

We will also demand that this construction solves the hierarchy problem, i.e. r_ ^ 37/, 
and thus for fixed ratio f 1/^2 the parameter m'^ (and hence A;_) should be accordingly chosen 
while we treat ^_ (or r_) as a free parameter. In principle C2 can always be chosen such 
that one obtains for ^_ an arbitrary value in the interval [1, 00] (with Aqo > 0) but we will 
assume ^_ to be smaller than fi/f2 in order to obtain a negative A;_ 

fc- = -logf^^) . (53) 

About the approximations we employed, notice that we used ( IT3l) as scalar potential and 
nowhere the assumption 0^ ^ M^ that needs to be fulfilled in order to make contact with 
the superpotential method. The sole constraint for the applicability of our results so far is 
flTSj) . This leads to the inequalities [using fl20|) ] 

m^/V^ < m^Pvf ^ -Alk^vl < 2AM^ , (54) 

and 

Z20'2 _ ^2 ij^^_^ ^ if^^^^ _ ^))2 ^ ^2 ^if^_^_ ^ ^^^^^ _ ^_))2 ^ 24M3 . (55) 

Depending on the parameters these constraints can be much weaker than those employed in 
the literature to ensure smallness of backreactions (g, [lO(] namely 0^ <^ M^. However, they 
still guarantee that the deformation of the CFT induced by the operators corresponding to 
the bulk scalar can be treated perturbatively at all scales up to the 4D Planck scale ^. 



^ In the case A;_ > these operators become strong in the IR which is potentially more cumbersome [10[ • 
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Notice that even if the conditions (IMI) and (155|) are fulfilled (which ensures that the 
relative error in A' is small) there still can be sizeable (cumulative) deviations in e^ since 
y4(ro) ~ 37. Nevertheless this will not affect the equations of motion of (p that only depend 
on A' and not on the exponential warp factor. Hence the expansion parameter deduced from 
the equations ( H3l) - (H^ is reliable. Still the detuning parameters need in principle (small) 
corrections to reproduce the correct hierarchy. In the following we call the regime where 
our constraints are saturated and our approximation becomes unreliable the regime of large 
backreactions. 

Let us compare these findings with the results from the superpotential method and those 
in the GW framework. The GW potential can be written as 

Vaw = r^vli—n [{A-lk_){l-i)^ + lk_]+const , (56) 



and the extrema fulfill 

4 2(4-/fc_)(l-e±) 



lk.^± (4-/fc_)(l-e±)2 + /fc_ ■ 
This leads to a quadratic equation for C,± with solutions that are of the form 



(57) 



for an appropriate C2. This coincides with the positions of the extrema in the detuned case 
found in (HSj) and fH9l) under the assumption fl50|) . The value of the potential in the extrema 
is given by 

VGw{r±) - Vgw{oo) = 2vll-\-^''^'\lk.)U{l - U) , (59) 

which is proportional to the difference in expansion parameters {K± — Aqo). The value of the 
radion potential can indeed be inferred from the expansion parameter A as we will see in 
section HV CI In conclusion our results agree [at least to order (9(0^)] with the superpotential 
method for the choice 

^_ = 1 (superpotential method) , (60) 

while the GW potential (up to a constant piece) is reproduced for the choice 

^_ ~ 1 + |/A;_|i/V2 - lk-/4: (Goldberger-Wise) , (61) 

under the assumption ( 150|) . Hence our approach unifies the results obtained both in the GW 
and superpotential approaches. 
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Finally some comments on the role of the original GW potential are in order. In our 
framework the potential obtained by Goldberger and Wise merely corresponds to a specific 
choice for the parameter ^_. From our point of view there is no special significance to 
this choice and it just results from imposing initially the same tension on the two branes 
when the backreactions are ignored. Once the behavior of the scalar field is determined 
one would be forced to adjust these brane tensions in order to cancel the contributions 



from the backreactions anc 
already mentioned in ref. 



to obtain a vanishing expansion parameter at the minimum (as 
6|). Hence this parameter choice is in no way a distinguished 
one. In particular the fact that the potential difference 1^^ scales with [IkJ)^/"^ results from 
the peculiar choice in (16T1) (or equivalently for the brane potentials). This scaling is also 
reflected in the radion mass as we next discuss. 

B. The radion mass 

In this subsection we will determine the radion mass in the detuned case along the lines 



of ref. ^5| while the decoupling of the linearized Einstein equations is demonstrated in 
appendix |Al It turns out that in the case of an expanding Universe the impact of the 
expansion parameter A can be absorbed in the radion mass parameter 

Kad = "^lad + 6A . (62) 

However the contribution from A is anyway negligible (being smaller by a factor e^^) leading 
to the equation [see eq. ( 1A15I) 



F" + 2A"F - 2NF' - V^F' = -ml,a e"'^F , (63) 

00 

which can be solved to obtain rrirad- If backreaction is neglected, A" = 0, this allows for 
a solution with F = 1 and m'^ad ~ 0- -^^^ small backreaction, F = 1 + f, this system of 
equations can be linearized as 



f" -2A'f' -2-^f' = -2A" 



m. 



rad 



'0 

2k^ 



g Jl-ml,,e-^^ , (64) 

while the boundary conditions in the limit 7^ — )■ C)0 read /'(O) = j'ij"^ = 0. The main 
difference with respect to the standard solution presented in ref. (3|] comes from the factor 
0q/0o that deviates from /c_ close to the TeV brane. 
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The solution reads 

fir) = 0'^(r) e^-^W £ dx (^^e-^^(^) - ^"~'^) ' ^^^^ 

that by construction fulfills the boundary condition /'(O) = 0. At this point we will neglect 
backreactions in A that would lead to corrections of order 0(0^) on the radion mass ^. Then 
the constraint /'(^o) = determines the mass by the equation 

For the scalar field solution in the limit of small backreactions flT^ . which can be written 



/,'(r) = vik.e'"-' (1 - g £('=+-'=-)('-"»)) , q = k+/k^ M - -\ 



(67) 



the above integral becomes (using lk+ = 4 — Ik] 



^ix/l i^{A-2lk-)ro/l 



ro 

For g > 1 the integrand actually contains a pole. However a more sophisticated analysis 
shows that the solution (!65|) is still regular due to the pre-factor 0'^(r) and that the naive 
integration is justified. Therefore one concludes that the radion mass is 



f-^ 



Observe that in the case without detuning (corresponding to a scenario that can be 
treated with the superpotential method) 

«- 1!^ • (™) 

the masses turn out to be 

2 f^ 

mrad,± ~ T^(^^-)' e"''*^' J^ (superpotential) , (71) 



When approaching the region of large backreactions one would have to take the above mentioned cu- 
mulative effect in the warp factor into account and the final radion mass would involve the corrected 
warp factor instead of the plain one. Notice that the determination of the radion mass only serves as an 
additional check on our radion effective action and the tunneling analysis of the radion is generally not 
affected by this subtlety. 
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which agrees with the physical radion mass fn^ad found in 15| and also with the interpretation 
that C,- denotes a minimum of the potential while C,+ is a maximum. Moreover in the case 
of GW, see eq. (l6Tll . the mass is 

<ad^^\lk-f/'e-'^'>/'^ (Goldberger-Wise) , (72) 

which indeed scales with |//i;_p/^ as expected from the GW potential in eq. fll3p . 

Notice that the radion mass is slightly below the TeV scale what is essential for an effective 
four-dimensional description. To compare the radion mass obtained here with the one from 
the radion potential involves the kinetic term which is the topic of the next subsection. 

C. Kinetic term 

So far we have only discussed the occurring expansion parameter A for the different 
solutions of the equations of motion. The main motivation was the relation ( IT^ that implies 
that the action is proportional to the expansion parameter. In the current subsection we will 
make contact between the five- dimensional system and the effective action of the radion. 



Let us start with the kinetic term of the radion. 
and here we briefly review the calculation of ref. 



t was derived in several ways 



ay 



211 



15|. We will use the metric ( JAip but 



will neglect the contributions from the scalar field. These effects are suppressed by a factor 



[Ik^Yvf/M^ 151] which is small under the constraints ( l55l) . The Einstein equations are then 
solved by the radion Ansatz F{x, r) = e^'^^^Ri^x) and the Einstein- Hilbert action contains a 
kinetic term of the form 

Krad = 6MH j d''x^\det-g,,\ (e^'-o/' - l) {dRf . (73) 

The correct normalization is obtained by the observation that in this background the geodesic 
distance of the branes is given by 

pro 

r{x) = dr{l- 2e^''^R{x)) = ro - I e^''/^R{x) , (74) 

Jo 

from which {dRY = l~'^e~'^'^°^\dry ~ l~'^e~^'^^\drY. Hence Krad becomes 

Krad = 12MH-' fd'x^f\d^t^\^{drfe-^'-/' 

= 12{Mlf J d^x^\detg^,\ ^ {d^if , (75) 
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with the notation fi = l~^e~^/K We do not absorb the pre-factor 12(ML)^ in the definition 
of /i because it will prove to be useful in the analysis of the phase transition ^. 

In ref. [g] the action resulting from the bulk scalar was right away interpreted as the 
negative potential seen by the radion. In the present case this could hardly be true since 
the action is positive for an expanding Universe. We better expect this to correspond to a 
maximum and not to a minimum of the potential, since on the other hand the radion mass 
is negative in this situation. Moreover an expanding Universe without any bulk scalar would 
also have a non-vanishing action, even though the radion does not see any potential in this 
case. 

In order to obtain the potential seen by the radion we must separate the action into the 
contributions from the expansion and that from the radion interacting with the bulk scalar. 
The contribution to the action from the expansion is according to (JTj) given by 

^exp = 12M^ f rf^V|det(7Miv|e-'^A , (76) 

and, considering the action provided by ( TT2l) . it results in a radion action 

Srad = -QM^ J d'x^\detgMN\e-^^A , (77) 

which must arise from the value of the radion potential at its extrema. 

As mentioned before an additional problem arises from the fact that the solution is not 
an extremum of the action in general relativity. To overcome this problem one can confine 



the system to a box, t G [—T, T] , and add the so-called Gibbons-Hawking 18|] term to the 
action ^^ 



Sgh = 2K , (78) 

JdM 

where K denotes the extrinsic curvature and dM is the boundary of the space-time manifold. 



Using the metric given in ([2]) one can evaluate the Gibbons-Hawking term [22 1 

Sgh = -18M3 J d'x^\detgMN\e-'^A , (79) 



^ Notice that this kinetic term was obtained by decoupUng the Einstein equations and the result is a factor 
of 2 smaller than the one from the more naive approach derived in [y, llOf and subsequently used in the 
analysis of the phase transition in lOJ, |ll|, |l3[ . 



^° Inclusion of this term is actually not relevant since the radion potential could be inferred from (frf]). but 
it will elucidate the comparison with ref. '^ later on. 
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and the total action is given by 

Stot = ^exp + Srad + ScH = -12M' f rf^X V^d^t^^C'^^A . (80) 

After integration over the fifth dimension the effective action for the radion hence reads ^^ 
Seff = 12{MlY f d^x^\detg^,\(^{dfi)'' -Vradil^)] 

+Ml I rf^x-^|det^^,|/?(4) + Sgh , (81) 

where we defined the four-dimensional Planck mass, {MplY = (Ml)^, and identified the 
radion potential in its extrema as 

Vraci{r±) = l-^A±/2 . (82) 

This is the essential relation that we use to connect the value of the radion potential to the 
expansion around its extrema. 

Let us compare this result with the one obtained in ref . 3| . There an effective action was 
obtained by perturbing the background metric but without taking a stabilizing mechanism 
into account. Their result reads in our notation (a is the scale factor) 

dt a' (h'e-'^'/H-' - ^^) ' - Vradir)] . (83) 

The last equality is obtained by partial integration which cancels the Gibbons- Hawking 
term at a time-like boundary as discussed above. The radion potential Vradi'f') was not 
specified in [3] but comparison with our result fIST]) confirms that the radion potential and 
the expansion (including the Gibbons- Hawking term) contribute equally to the action in the 
extrema of the radion potential. 

D. Interpolating potential 

In the following we will approximate the effective potential seen by the radion. The 
equations of motion have only stationary (up to the Hubble expansion) solutions in the 



^^ Note that we integrate over the orbifold S1/Z2 and hence twice over the bulk as it is customary in the 
hterature on RS models. 
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three extremal situations r G {r_,r_|_,oo} [or ^ G {^-,^+,0}] and we use these three values 
to provide a fit to the potential. The positions of the extrema of the radion potential are 
given by eqs. (H5ll and (HHl) . We parametrize the potential by 

KUOocr-'(e-e-)(e-e+), (84) 

and integration yields 

VraAO OC r (^ - ^^^4t^ + ^-^) ' (85) 

\u + 2 w + 1 i^ / 

where the integration constant is fixed to zero by K-ad(O) = 0. The parameters u and the 
pre-factor can be adjusted in order to reproduce Vrad{^±) given by eqs. ( 15T|) . (l52l) and ( l82ll . 
This corresponds to the expression 

Kad(0^' = Ae-^+'-p(0, (86) 

with 

24 M3 + ' ^ ^ 

and 

„ 1 - ?- / ? A "+'*"^ ^(^ + 1);' - u(u + 2)(;_ + ^^)^ + (c^ + i)(c^ + 2);,^^. 
"^-^"Irj -K--W-2«, ■ '**'' 

The function P{C) is normalized to P{^-) = (.^- — 1) and u is determined by P{^+) = (.^+ — 1), 
which yields 



.^-7 u;(e--e+)-2e+ 1-e-' ^^^^ 

and hence u ~ —k^/k_. 

Another ingredient useful for the fit is the radion mass at the extrema where the second 

derivative of Vrad turns out to be 

Using the standard kinetic term for the /i field 

^ = /-I e""/' , d^/dn = -Ik^^/fi , (91) 

one gets the expression for the mass 



m^ 



KM^ (1)^ 



14- - 4+J ± 24±/w 
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If ^_ is not too close to unity the last factor becomes {^± — 1) in agreement with (169|) . On 
the other hand in the superpotential limit fl39l) 

e-^1, e+^l-2/a;, (93) 

the last factor behaves as [using (!89|) ] 



in agreement with eq. fl7T]) . This shows that our fit is indeed consistent with the effective 
action (IHTI) . 

The above parameterization is reasonable as long as ^+ < 1 which is the case we will 
consider hereafter. For the ^+ > 1 case on the one hand eqs. (1^ and (1S21) imply A+ < Aqo, 
so that Vrad{^+) becomes negative and at the same time transforms into a minimum according 
to (l69l) . In this case our Ansatz for the metric ([2]) does not allow for additional solutions 
between the two minima at C,± while according to (!69|) any additional solution should be a 
minimum. Therefore our Ansatz in ([2]) produces in this case two local minima of the radion 
potential in configuration space but no local maximum. This complicates the question of 
what the radion potential might look like in the case ^_|_ > 1 so that we will disregard this 
case in the following. 

V. HOLOGRAPHIC PHASE TRANSITION 

In this section we will present the discussion of the holographic phase transition at finite 
temperature along the lines of refs. 10|, llll, ll3| ^^. At finite temperature the system allows 
for an additional gravitational solution with a black hole singularity in the bulk. This 
AdS-S metric describes in the AdS/CFT correspondence the high temperature phase of the 



system 2J, |25|]. This phase starts dominating at temperatures of the order of the TeV scale. 
In fact the potential difference between AdS-S and pure AdS phases is given by [10|] 

-An^Mim^^^^ , (95) 

where T^ is the scalar field parameterizing the distance between the horizon and UV brane 
and T is the temperature of the system. On the other hand the difference between the 



12 



A precursor at zero temperature can be found in |23j . 
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RS-GW and pure AdS phases is expressed by eq. (159|) which should equal eq. ( p5|) at the 
critical temperature T^ of the phase transition. It turns out that typically T^ is between the 
electroweak and TeV scales ^^. 

Starting from a hot universe in the AdS-S phase a first-order phase transition towards 
RS-GW may happen below Tc . In the five-dimensional picture this means that the IR brane 
emerges from the black hole horizon and to minimize the bounce action this crossing has to 
happen far from the UV brane. For this reason it is commonly assumed that the bounce path 
consists in moving the horizon away from the UV brane till arriving to the (unstable) pure 
AdS phase, and subsequently displacing the IR brane from r = cxd to r = r_. This reasonable 
assumption fixes the bounce path and it reduces the study of the tunneling probability to 
the usual analysis of the one-dimensional bounce J26|, |27| in which the bouncing scalar field 
is identified with /i (T/j) in the part of the path between pure AdS and RS-GW (AdS-S) ^^. 

In particular the path in the AdS-S space can be simplified by observing that the kinetic 
term of the field Th has a small pre- factor ^^. Consequently this part of the path becomes 
extremely short once Th is canonically normalized, so that in the AdS-S region the potential 
seen by the bouncing field can be approximated by a step function. 

We will release the radion field fi from a certain initial position /iq and we will evolve it 
to the point ft = (corresponding to the pure AdS phase) according to the 0(3) bounce 



equation [26|, |27| 



a> + 2^ = d.Vra, , (96) 

where p^ = x'^. In the bounce solution the radion field should arrive at p = with the 
kinetic energy necessary to jump and stop on the top of the AdS-S minimum 

4vr^T^ = 6(9,/i)^|^^„ , (97) 

and this solution is used to determine the bounce action S^/T. 



^^ In reporting the expression (|95p we omit the subdominant contribution of the bulk matter fields since this 
correction plays a minimal role in our analysis. The same consideration holds for the thermal corrections to 
the RS-GW potential. We assume nearly all SM fields to be fundamental so that their thermal contribution 
to the two-phases free energies is similar and thus it cancels out. 

^■^ For more details about the AdS-S solution and the bounce of this phase transition see refs. 1CH13|. 



^^ The estimate of this factor is controversial. We checked in the numerical evaluation that admitting a 
sizeable kinetic term would only lead to slightly smaller tunneling temperatures. 
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In order to explicitly calculate the tunneling probability it is useful to rewrite the radion 
action flHTl) as 

Seff = 12 (Mlf I d^x^/\deigMN\ T^A' - VradilA , (98) 

where the potential, neglecting a constant piece, is 



4 u{u + 1){UJ + 2) [ fl 



u 



(l-li)-2|^VA 



X 



1 ^M-"- 1 f^Y'-(, + L]^LL 



with u determined by eq. ( 189|) and 



(99) 



j-l I ^2, 



-i/(/fc_) 
/i_ = 1-' (^^e-J , (100) 

A = ^^(/A:-)(/M(^-l)e^- (101) 

The dimensionless tunneling action will only depend on the parameters A, u and fi^/fi^ 
and not explicitly on /i_. This can be seen by using the conformal transformation 

x^ -^ a^^x^, Vradil^) -^ a'^Vradiafi) , (102) 

which is equivalent to a rescaling of all dimensionful quantities, in particular fi_ — )■ fi_/a. 
Therefore the relevant scale involved in the bounce solution is /iq and the functional de- 



terminant 



26[ of the tunneling process has to be proportional to Hq. Tunneling can hence 



occur for bubble action Sb 

St ~ log ^ ^ 4^ + 41og ^ < 140 , (103) 

where it is used that the expansion parameter is the cosmological constant of the RS-GW 
phase (l52l) . Observe that if the tunneling takes place during radiation domination the Hub- 
ble expansion is quadratic in temperature and the right-hand side of this equation would 
increase when the tunneling temperature decreases. However in the present case the Hubble 
parameter is dominated by the vacuum energy (which can even lead to a short period of 
inflation [l3|, |27|) and a smaller temperature in fact lowers the right-hand side via its depen- 
dence in the functional determinant. Notice that even though the four- dimensional gravity 



25 



contributions to the action are significant (due to the expansion) gravitational effects ^^ in 
the tunneling process should not relevantly modify our results 28 |. 

As long as the release point is far away from the minimum and maximum of the potential 

(^^/^_)i/l^-l « ifi « 1 , (104) 

the field experiences a nearly conformal potential. If ^_ is not too close to unity the position 
of the maximum, given by 

(«./«-)-"- < (^«-)"'" = e-'-" = I,- = #^ , (105) 

\V2 J Jr^lanck scale 

is such that there is a large hierarchy between the position of the minimum and the maximum 
and for a large range of release points the potential is nearly conformal. 
In the conformal case the potential is of the form 

Vconf -^ -/t/i^ • (106) 

Solving eqs. ( l96l) and (1971) with this potential gives an action, temperature and bubble size 
for the 0(3) symmetric bubble 

S^/T ~ 217.0 K"^/^ (M/)^ T//io ^ 0.103 fi:^/^ p/^o ^ 3.45 R"^''^ . (107) 

However the potential (I106p is normalized to the origin instead of being normalized at the 



AdS-S minimum 26| and we have to add to the action the omitted contribution 

levr^ 



T^p^ {Mlf ~ 72.3 R-^'^ {Mlf . (108) 



3 
This yields for the tunnel action 

S^/T ~ 289.3 R-^'^ {Mlf . (109) 

Besides thermal fluctuations the potential barrier can also be overcome by quantum 
fluctuations. In this case the bounce solution is 0(4) symmetric and eq. (p6l) has to be 
replaced by 

5> + 3^ = 9,Kad. (110) 



^^ Gravitational effects become important when nucleated bubbles are of order l/vA- 
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This tunnel configuration is relevant for temperatures that are below the inverse bubble 



radius 
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271]. In the present model the bounce solution results to be in this regime and the 
quantum tunneling competes with thermal tunneling. In the conformal case, the solution to 
the bounce equation is given by 



2/^0 
2 + Kp2/io 



M = TT^^ > (111) 



which leads to the action 



27r2 
^4 = 12{Mlf—- , (112) 



while the corresponding temperature vanishes. Small deviations from the conformal case will 
only lead to a relatively small temperature while sizeable temperatures can only be obtained 
when the conformal symmetry is broken by a release point which is not too far away from 
the minimum of the potential. In this regime our nearly conformal approximation (1112p 
underestimates the real action that could be reliably obtained in the thin- wall approximation. 
Still one can generally conclude that for the same release point the thermal tunneling leads 
to larger temperatures than the quantum tunneling since the friction term in the bounce 
equation is smaller. Hence if both tunneling modes are feasible the system tends to tunnel 
by thermal fluctuations. 

In the near conformal case a good approximation is given by using a conformal potential 
normalized at the release point 

Vnear-conf = Vrad{f^o) I 1 , (113) 

such that 

_^_Ka^_ (114) 

In this approximation the minimum of the bounce action as a function of fiQ (or equivalently 
T) is given by the maximum of Kad(yUo)//^o- According to fl85|) the function VradifJ') / fJ'^ is 
nearly polynomial in terms of ^ with one extremum between ^_ and ^+. The derivative of 
this function at the endpoints is given by 

4 Vradifi) 



d^ = d/-'^^^ 



^^^ 



A reasonable estimate of the position of the minimum is 



= ^^(^M(e±-i). (115) 



d.-d. ^e-+e.-l, (116) 
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with the approximate value 



'd^-d^ 2 



12 M3 



(/M(e--l)(l-e4 



;il7) 



For example a nearly maximal value is given for ^_ k. 1.5, t'i/f2 ~ 3. In this case the 
coefficient is given by k ^ 9.2 x l{)^^v\/M^. For this choice of parameters and with vi = 
4M^/^, A^ = 3, k^l ~ —0.019, which create the correct hierarchy, a comparison between this 
approximation and the full numerical result is shown in fig. |2l It turns out that the estimate 
is quite reliable and in particular no appreciable difference is found for 5*4. 

The bulk cosmological constant is related to the five- dimensional Planck mass by the flux 
A^ of the background fields which, according to the AdS/CFT correspondence, translates 
into the rank of the gauge group of the CFT in four dimensions. In the present context this 
relation is only known to leading order in A^ and here we use (Ml)^ = {N"^ — l)/167r^ as 
definition 29|]. Then considering the maximal allowed h for the thermal tunneling action 
one obtains the bound 

-3/4 

(118) 



Ss/T > 61.5 X (A^^ - 1) ( J^) 
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FIG. 2: The left plot shows the comparison of S's/T and 5*4 as a function of the release point ^uq 
between the approximation in (J113p and the full numerical results. The two curves for 5*4 lie on top 
of each other and significant deviations only occur for a release point very close to the minimum of 
the potential. The right plot shows the temperature T as a function of the release point /io- The 
used values are vx = AM^^"^, v^ = vi/3, iV = 3, ^_ = 1.5 and kA ~ -0.019. 
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and a nucleation temperature that is substantially smaller than the scale //_. Using the 

criterion f ll03p for tunneling this translates into 

/ 2 \ 3/4 
iV' - 1 < 2.3 X (-^j (thermal) , (119) 



which is quite similar to the bound obtained in 10|. Analogously for the case of quantum 
tunneling one finds 

^4 > 54.3 X (7V2 _ 1) (^<^ , (120) 

and 



iV^ - 1 < 2.6 X ij^j (quantum) . (121) 

Depending on the value of vf/M^ the weakest bound between (I119p and (I12ip constitutes 
the minimal requirement on A^ allowing for a first-order transition. 

Alternatively one can saturate the bound of applicability of our analysis (155|) which reads, 
for ,^_ not too close to unity, as 

4f^--l-^V<4, (122) 



M3 V k+J ^ 2P 

such that 

R < 1 (^--i)(i-^+) < i n23) 

'^~2(e--l-/t_/A:+)2~2' ^^^'^' 

which gives for thermal tunneling 

A^ < 6.8 (thermal) , (124) 

while for quantum tunneling one finds 

A^ < 12.5 (quantum) . (125) 

Thus there is a small window where quantum tunneling is possible while thermal tunneling 
is not. 

We conclude that even though the limits on tunneling derived including backreactions 
are parametrically similar to the results obtained without them our anal ysis shows that the 
constraint on the parameter space used in the literature, 0^ <^ M^ [6|, Il0l-ll3|. is actually 
more conservative than necessary. 
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VI. GRAVITATIONAL WAVE OBSERVATIONS 

The main input parameter for the determination of the gravitational wave spectrum is 
the inverse duration of the phase transition that, normahzed to the Hubble parameter, is 
given by 

This equation assumes that the temperature is proportional to the inverse scale factor. 



T oc exp^—y/Xt), which is also true in the AdS-S phase 29|. In this section we will present 
analytic as well as numerical estimates for /3/a/A. 

The dependence of the coefficient R in the near-conformal case is approximately polyno- 
mial at the release point ^o- Besides the tunnel action should grow to infinity if ^o approaches 
either ,^_ or ^_|_ ^ 2 — ,^_ . A reasonable fit is then given by 

Ri^o) = ^^ma J I - (^^) J . (127) 

Hereby and depending on the model parameters the value of ii^ax is somewhere in the range 

3.1 X 10-^{N^ - 1)4/3 < R„,ax < 0.56 (thermal) , 

3.6 X 10-^{N^ - 1) < i^max < 0.56 (quantum) , (128) 

where the lower bound arises from the constraint that the system will tunnel [eqs. (I109p 
and f lll2p ] while the upper bound results from a numerical equivalent of (11231) obtained by 
exploring the parameter space constrained by the requirement of small backreaction ( l55l) . 

For the thermal tunneling action fllOOp the inverse duration of the phase transition is 
given by (notice that approximately /iq oc T as shown in fig. [2]) 

= |fH«:-)^%^7l-^. (129) 

At the first sight it seems that /3/a/A can become arbitrarily large for ,^_ — ?■ 1. However 
if ,^_ approaches the unity the constraint (!54l) becomes more severe than (!55l) and Kmax is 
according to (11171) bounded by 

^ ''':(^M(^-l)(l-e+)<7^4(^-l)(l-^+)' (130) 



12 M3' +^VS ^V -^^^ lk_y2 
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such that the maximum possible value of 13 /yK in fact becomes constant in this limit. Hence 
one can obtain a conservative estimate by assuming that ^_ is in the transition region where 
the two constraints (|5^ and fl55|l are equally severe 

e- - 1 ^ ^nk:\- , (131) 

and hence 



Sz V2 Umax L >^ 



P/VA<3y/-lk_^—-^^^Jl-- . (132) 



Noting that 



one obtains the bound 



v/=^- ^ y^^lkle'-^- < J-L- , (133) 

Vi V 2er_ 



P/VA <xI^ SsRrna^ I^ ^ 30^^./l--^ . (134) 

Z 6 7 _ 1 ri \ i^m.ax f^ V f^m.n.x 



Eq. (11281) for the thermal case implies 



f^max ^ loU 



< 



K ^ {m- 1)4/3 

which finally yields a bound on /3/vA. 

Analogously for quantum tunneling one finds 



(136) 



^/VA<,/lis,5=../l-^(^^). (136) 

The last factor is typically smaller than in thermal tunneling (where it is unity, cf. fig. [2]) so 
that commonly quantum tunneling happens only when the system cannot tunnel by thermal 
fluctuations, which implies 

/i4n \"^/^ 

Rmax < (i^ievr^J {N' - ly/' . (137) 

Besides the largest R that can be realized in the present model is generally given by Umax < 
0.56 according to eq. (I128p . On the other hand the quantum tunneling condition (I128p 
implies 

K > (2 X 140)"^ (iV^ - 1) , (138) 

such that successful quantum transitions which are not spoiled by earlier thermal tunneling 
lead to the constraint 

^^ < min [0.85 (iV^ - 1)^3^ 155 (jys _ i)-i] . (139) 
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The first bound increases with A^ and arises from the requirement of no thermal tunneling 
while the second bound stems from viable quantum tunneling and decreases with A^. As a 
result, as far as A^ decreases the quantum tunneling bounds on /3/vA become first weaker 
(while thermal tunneling becomes less likely) and then stronger (as thermal tunneling is in 
general not possible for N > 7). Finally for A^ > 12 tunneling is impossible altogether. 

The approximate values and numerical results on the bound on /3/a/A are given in Table [H 
Already for N = 2 one observes a rather strong phase transition (a > 1). We see that for 
small A^, when the release point is close to the minimum of the potential, our approximation 
overestimates the numerical result for /3/-\/A, while for larger values of A^ it yields fairly 
precise results. 

The second ingredient for the gravity wave spectrum is the vacuum energy normalized to 
the radiation energy of the system (traditionally denoted by a in the literature) 

12(M/)3(V;,,(oo) - Vradir-)) l-'A 



a 



(140) 



Pradiation -^^ -^ n 

In the present system this parameter is much larger than unity as it is shown in Table [T] due 
to the large supercooling. This also implies that the nucleated bubbles expand with near 





thermal 


quantum 


N 


Papprox / V A 


/5n«m/\/A 


Oinum 


l^num/VA- 


Oinum 


2 


< 1230 


<770 


>3.2 


< 15 


> 1.0 X 10^^ 


3 


<235 


<315 


> 10 


<33 


> 5 X 10^ 


4 


< 131 


< 143 


>50 


<45 


> 2.5 X 10^ 


5 


<62 


<67 


>800 


<56 


> 4 X 10^ 


6 


<29 


<30 


>105 


<63 


> 1.3 X 10^ 


7 


<6.5 


<6.0 


>108 


<71 


> 5 X 10^ 


8 




<54 


> 5 X 10^ 


9 




<40 


> 8 X 10^ 


10 




<27 


> 3 X 10^ 


11 




< 17 


> 4 X 10^° 


12 




<5.4 


> 5 X 10^4 



TABLE I: Upper limits on /3/v A and lower limits on a for all possible values of N. 
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luminal velocities. The vacuum energy is very efficiently transformed into bulk motion of 
the plasma or directly into kinetic energy of the bubble wall and resides in a very thin shell 
around the bubble wall [30|. In this sense the phase transition is extremely strong. 

The main mechanisms of gravity wave production during a first-order phase transition are 



bubble collisions 
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ll-l36|. turbulence 37H40l| and magneto-hydrodynamic turbulence 4ll. |42|. 
In our analysis we will focus on the production mechanism by bubble collisions for the 
following reasons. First, while reliable information about the gravity wave spectrum in 
the case of bubble collisions is provided by computer simulations, the analysis of turbulence 
relies typically on additional assumptions {e.g. on the overall normalization of the spectrum). 
Second, the peak frequency of the spectrum produced by turbulence compared to the one 



by bubble col 
fiuid motion 



isions is typically suppressed by the eddy velocity that occurs in the turbulent 
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38|. At the same time the spectrum for turbulence falls off faster than 



the one by bubble collisions. Hence even if the different components are of similar size they 
might be disentangled because they lead to a double peak structure or at least a knee in 
the spectrum. Finally, recent studies have shown that very strong phase transitions most 
probably have a runway behavior of the wall |43[ and in this case most of the vacuum energy 
is transformed into gradient /kinetic energy of the Higgs and not into collective bulk motion 
of the plasma 30|. This should also reduce the portion of energy that leads to turbulent 
plasma motion. In the following we will focus on the spectrum produced by bubble collisions. 
In the above discussed limit, a ^ 1, the energy fraction for colliding bubbles of the 

2 



2ri* 



h'n 



7.7 X 10" 



A 



(3 



(141) 



In a standard cosmology this energy density is diluted to today's observed energy fraction 

2 



h^Q = 1.3 X 10" 



/3 



(142) 



The only assumptions that enters here is that the Universe is dominated by radiation after 
the phase transition. Immediately after the phase transition the energy fraction stays con- 
stant up to the time of matter-radiation equality when the energy fraction of gravitational 
radiation starts to become suppressed together with all the other light components of the 
plasma. 
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The frequency peak of the spectrum at time of generation is in the hmit of very strong 
phase transition, a :» 1, and large wall velocities, f^ ~ 1, given by [36 1 

/* = 0.23/3. (143) 

The red-shift to today's observed spectrum depends on the reheating temperature after the 
phase transition according to 



/ = 0.23/3 



T, 



reh 



(144) 



where Tq is t 



is of form 



le observed temperature of the cosmic microwave background. The spectrum 



36| 



n{f) = Q: 



3.8 r-«/ 



2.8^1.0 



(145) 



1.0/3-8 + 2.8/3-8 

Several sample spectra are given in Fig. [31 

As long as the expansion is negligible during the phase transition the reheating tempera- 
ture Treh can be determined by energy conservation. Since the phase transition is extremely 
supercooled at the typical nucleation temperatures the energy density of AdS-S phase is 
comparable to the one of pure AdS and thus 



n 



9*^Veh ^ l2{Miy{VU0) - V;.rf(/i-)) = 6 (M/)3/-^A 



(146) 




10000 



FIG. 3: Several example spectra of gravitational waves. The straight (dashed) lines are for a 
reheating temperature Tj-e^/vA = lO"*^^ {T^-eh/v^ = 10"*^^). From bottom to top the plots use 
/3/\/A = 11000,300,100,30,10). The sensitivities of the LISA and BBO experiments are taken 

n 

from 14411. 
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where A = Aqo and (?* denotes the effective number of degrees of freedom just after the phase 
transition. Notice that for g^, < 1207r^(M/)^ the critical temperature, given by 

A-K^T^ = 61-^ A , (147) 

is smaller than the reheating temperature and percolation is followed by a period of phase 
coexistence ^^. Nevertheless we will assume that the red-shift of the gravity waves is given 
by the naive expression involving the reheating temperature fll46p . 

We now give a conservative estimate for the range of possible values for Treh/vA in a 
realistic model. According to fl52|) and ( 169|) the expansion parameter is related to the radion 
mass by 

A^^m^e-^^-/', (148) 

such that together with the definition of the four- dimensional Planck mass, {MplY = {MlY, 
one obtains 

9.^^ = 484^e-/^ . (149) 

Introducing the TeV scale fx- = l~^e~^~^^ and the relation {Mpiy = (N"^ — l)/167r^ the peak 
frequency can be written as 

/ = 1.77 X 10-3 mHz -^(N' - 1) 1/4/^1^ . (150) 

VA^ ^ V TeV^ ^ ^ 



Realistically the last factor varies within a range of [0.3, 10]. The corresponding ranges for 
the possible positions of the peak structure of the gravitational wave spectrum compared 
with sensitivities of future experiments are shown in Fig. HJ For A^ < 3 the phase transition 
could be weak enough to prohibit detection by LISA. Besides for these values the 1/A^ 
expansion is clearly questionable and stringy loop contributions invalidate our analysis. On 
the other hand for A^ > 4 the model would lead to a gravitational wave signal that is 
observable by LISA. 

VII. CONCLUSION 

Motivated by the sizeable impact of backreactions on the radion mass determined with 
the so-called superpotential method [l5|, [l6| , we studied backreactions to the radion poten- 
tial in a perturbative scheme. Our approach includes backreactions on the five- dimensional 



^"^ Assuming that some SM fields be composite relaxes the condition on g^. 
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f/mHz 



FIG. 4: The regions denote the possible positions for the peak of the gravity wave spectrum 
depending on the parameter N. The signal will be detected by LISA, BBO or BBO correlated 
when it stands above their respective lines. 



gravitational fields from the bulk scalar under the constraints in ea s. (IMl) and (!55|) which 
are less severe than the usual assumption of weak fields, 0^ <^ M^ p, llCH13| . An immediate 
consequence of the backreactions is the necessity to generalize the metric to include the 
Hubble expansion of space-time in four dimensions. In our framework a judicious choice of 
brane tensions allows to reproduce the radion potential found by Goldberger and Wise as 
well as (within the applicability of our approach) the radion mass obtained with the super- 
potential method. We then clarify the apparent paradox that the radion mass computed in 
the GW and superpotential frameworks scale (for a bulk scalar mass m? ~ Ak_/l) as |/i;_|^/^ 
and |fc-P, respectively. On top of that in the regime of sizeable detuning of brane tensions 
we found that a scaling proportional to fc_ is possible. This results in a larger radion mass 
for fixed TeV scale and in turn it leads to a deeper radion potential which has significant 
impact on the thermal phase transition of the system. 

The phase transition constitutes a serious problem for the holographic interpretation of 
models with a radion stabilized by a bulk scalar 10l-|l4|: the tunnel action scales as A^^ with 
the fiux A^ of the background field such that the metastable symmetric phase becomes, for 
large values of A^, effectively stable and the conformal symmetry is never broken during the 
course of the Universe. One possible solution to this problem is to consider more general 
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scalar potentials (as e.g. done in ref. [45|) which modify the behavior of the bulk fields 
considerably in the IR. Analogously large backreactions can have potentially the same effect 
by decreasing the 'effective' N through a stronger warping close to the TeV brane. 

In our analysis of the phase transition, and compared to former work [lO|, lUl, we did not 



use the common thick or thin wall approximations to obtain the tunnel action but devised 
an approximation that is tailor-made for nearly-conformal potentials. Besides notice that 



we used the kinetic term obtained by decoupling the Einstein equations 15| which differs 



by a factor of 2 from the one used in 10|, lUl, ll3|| . The main consequence of our approach 



is that the tunneling action depends logarithmically on the temperature and can lead to 
a couple of e-folds of low scale inflation without tuning the model parameters (this was 



already observed in the numerical analysis of ref. [13|]). While a few e-folds of inflation 
cannot solve the horizon problem a low phase transition temperature has a large impact on 
gravitational wave production since, in this case, the energy stored in the vacuum bubbles 
during percolation at the end of the phase transition is many orders of magnitude larger 
than the energy stored in the thermal plasma. 

In summary we found that, in the regime of large backreactions, the deeper radion po- 
tential leads to a significantly weaker phase transition and numerically the absolute limit 
A^ < 13 applies. Besides we reanalyzed the gravitational wave spectrum produced by the 
first-order phase transition. We conclude that as long as stringy corrections can be neglected 
(specifically A^ > 3), the model leads to a stochastic background of gravitational radiation 
that can be observed by LISA (see fig. H]). 
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Appendix A: Linearized system of equations and the radion mass 



In order to calculate the radion mass we will follow the analysis of ref. 15| including the 



expansion of the Universe. For the perturbations of the metric we use the Ansatz 

(j){x,r) = (/)o(r) + (f{x,r) , 

ds^ = e'^+2^("''')^^,dx^dx'^ - (1 - 2F(x, r))^dr^ , (Al) 

where g denotes the induced (iS'4 metric as in eq. ([2]) 

g^^dxf'dx'' = df - e^^\dxl + dxl + dxl) , (A2) 

and solve the linearized Einstein equations, 5R^y = ST^j^y. The corresponding entries for the 
Riemann tensor are given by 



5R^.u -- 


= g^, i-DF + F" + WA'F 


SR,5 -- 


= -QA'd^F - 3d^F' , 


SR,, -- 


= -2nF - 16A'F' - 4F" , 



(A3) 
where we defined the d'Alembertian operator in curved space time {g = det g^u) 

UF = g-^'Mah^'^duF) = g^^d^d^F + sVAe^^^^tF . (A4) 

The energy-momentum tensor is given by 

5T,u = -^g,u{y\<i),)^ + 2V{<i),)F) 



■^9^^uY.\ 



^^'('^°)^ + 4A„(0o)FU(r-r.), 



(90 

5T^5 = 20o9^v5 , 

4 16 

5T,, = 40'o¥.' + -y'(0o)¥^-yn0o)i^ 



+ lY.{^-^^-2X.[^o)F)5[r-n)^ m 



and the scalar equation is given by 

ny,-y,"-4AV + ^(0o)¥^ 



- E (^^ - 2A.(0o)f) 5{r - n) - 60'oF' - A^-^F . (A6) 



and eqs. ©-(HO]) 

this can be brought into the form 
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The equation for R^^ can be integrated to yield 

ct>W = -\{F' + 2AF). (A7) 

Next consider the equation ^g^^^SR^u + ^-Rss in the bulk 

SDF + QAF' - 2AA'^F - 6A"F + 3F" = 8VF - 40'o(^' . (A8) 

With eqs. ([Ml) 

2 ./ ./ _ pn ^ 2^"F + 2A'F' -^{F' + 2A'F) , (A9) 

00 

- -V = A" + AA'^ - SAe-""^ , (AlO) 

DF = F" + 2A'F' + AA'F - 2-^ \F' + 2A'F] 

00 

+6Ae-2^F . (All) 

A minimally coupled scalar in four-dimensional curved space time fulfills the equation 

(e^^D + m2)F = , (A12) 

which leads to the equation 

(-m^ - 6A)e-2^F = F" + 2A' F' + 4A"F - 2^ [F' + 2A'F] . (A13) 

00 

Hence in an expanding Universe the sole difference on the equation for the radion mass 
with respect to the usual radion equation can be taken into account by a shift in the mass 
parameter 

m^ = m^ + 6A . (A14) 

Since the Hubble parameter is many orders smaller than the radion mass this shift can quite 
generally be neglected. Using the Ansatz F = e^^F this leads to 

F" + 2A"F - 2A'F' - 2^F' = -m^e'^^F , (A15) 

00 

which is the starting point of the analysis in section IIVBI 
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